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A Canonical Hamiltonian Derivation of Hawking Radiation 
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We present a derivation of Hawking radiation based on canonical quantization of a massless scalar 
field in the background of a Schwarzschild black hole using Lemaitre coordinates and show that in 
these coordinates the Hamiltonian of the massless field is time-dependent. This result exhibits the 
non-static nature of the problem and shows it is better to talk about the time dependence of physical 
quantities rather than the existence of a time-independent vacuum state for the massless field. We 
then demonstrate the existence of Hawking radiation and show that despite the fact that the flux 
looks thermal to an outside observer, the time evolution of the massless field is unitary. 
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Hawking radiation |l] is one of the most interesting 
phenomena encountered in pre- quantum gravity. It is a 
robust phenomenon in that it has been derived in differ- 
ent ways (see e.g. (2J |[ ||) and all derivations yield the 
same conclusion: a black hole of mass M emits (nearly) 
thermal radiation with a temperature Th = l/(8irGM). 
To the best of our knowledge however, no derivation 
of Hawking radiation discusses the problem within the 
framework of canonical quantization. This is one reason 
why the issue of whether the time evolution of the semi- 
classical theory is unitary has been a subject of debate. 
In this paper we present a simple canonical quantiza- 
tion procedure which leads to the usual results within 
the framework of a unitary quantum theory. Moreover, 
it enables us to compute the energy-momentum tensor 
as a function of position for all times, allowing us, as a 
matter of principle, to give a self-consistent discussion of 
the back reaction problem. 

The reason why it is not obvious that a canonical quan- 
tization scheme is possible in Schwarzschild coordinates 
is because surfaces of fixed Schwarzschild time are not 
globally space-like. In Lemaitre coordinates, however, 
surfaces of constant time are space-like and extend from 
r = to r = 00. In this coordinate system one can 
canonically quantize a field on a fixed time surface and 
then consider its time evolution from that point on. The 
immediate and most striking result of this approach is 
that the Hamiltionian of the massless field, computed 
in the background of the Schwarzschild black hole, is 
time-dependent. From this point of view, there are no 
uniquely defined eigenstates (in particular the vacuum- 
state) of the the quantum mechanical problem, in spite 
of the fact that we are dealing with a free field theory. It 
is this time dependence which is ultimately the origin of 
Hawking radiation. 

The metric of the eternal black hole has the following 



form in Schwarzschild coordinates: 



ds 2 = - 1 - - )dt 



1 ) dr" 

r - 
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where a = 2GM, G is Newton's constant and M is the 
mass of the black hole. Because the coefficient of dr 2 
changes sign at r — a surfaces of constant Schwarzschild 
time go from being space- like to time- like at this point. 
This does not occur in Lemaitre coordinates, A, 77, which 
are related to t, r by 
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A = t+2y/ra+a\n 



In these coordinates the metric becomes: 



ds 2 = -d\ 2 



-dr)' 2 + r 2 dfl 2 



(2) 



(3) 



and from this form we see that surfaces of constant A 
are everywhere space-like. Thus, these surfaces can be 
used to carry out the canonical quantization of the field 
theory. (Note, since r is a function of A, the metric is 
explicitly A or time-dependent and this translates into 
the time-dependence of the Hamiltonian.) 

The action of a minimally coupled massless scalar field 
in the background of the Schwarzschild black hole is 



S 



d^xy/^g^di^d^. 



(4) 



Since this action is rotationally invariant we can discuss 
one angular momentum mode at a time. From now on 
we will focus on the S-wave component of the field and 
define 4>o — 3>/\/47r. In Lemaitre coordinates standard 
manipulations lead to the Hamiltonian 
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where ttq is the canonical momentum conjugate to (j)o; 
i.e., 7r = a(j] — \)d<fio/d\. As always, 7r and </>o are 
assumed to satisfy the equal time commutation relation 
[7r o (A,?7),0o(A, V)] = -iS(r) - rf). 

The dynamics of the field </f>o and its canonical mo- 
mentum 7To are governed by the Hamilton equations of 
motion, which are equivalent to the statement 



d 



(V - A) 



dX 



d_ 

drj 



(r] — X)r 



0. (6) 



If we choose the initial space-like surface to correspond 
to A = 0, so that -q — 2/3 r^Jr/a and examine the A = 
Hamiltonian, we see that it is useful to rescale the field 
and its momentum as follows: 4>q = 4>i/r, ir = ^/rani. 
After this rescaling tti and <f>i satisfy the commutation 
relation: [iri(r), 4>i(r')] = —iS(r — r'), and the A = 
Hamiltonian takes the simple form 



n 



A=0 




(7) 



Clearly, this is just the Hamiltonian for the radial mode 
of a free scalar field and it is easily solved. As usual, self- 
adjointness requires that we impose the boundary condi- 
tion 4>i (r = Q) = 0, which implies that <f>\ and tt\ can be 
written as 



Although the assumption of large u> may seem unjus- 
tified, these two solutions do provide a very good ap- 
proximation to the exact solutions both at r/a — > 1 and 
r — > oo, which are the two regions that are of primary 
interest to us. (Note, however, that these solutions are 
invalid for r — > 0, a point we will have more to say about 
below.) Examination of these solutions shows that they 
are characterized by the fact that they are constant along 
ingoing or outgoing null geodesies. Abstracting this fact 
we define the geometric optics approximation as the as- 
sumption that the time evolution of the field 0o is given 
by 



<Po (A, r) 



l r~ 



MA + Si(r)) + fc(A + S a (r)) 



(10) 



where (f>\ and <fo are two arbitrary functions that can 
be determined from the requirement that at A = they 
can be written in terms of the field <po and its canonical 
momentum t:q. 

Imposing this requirement to determine the two arbi- 
trary functions ^^(Si^M) = fi,2( r ) we find 

M\r) = ^{f 1 (x 1 ) + f 2 (x 2 )}, (11) 



where 



Mr) 




sin(wr) (a+ + a u ) ; 



7Ti(r) = i I sin(wr) (a,. 1 , - n , 



(8) 



where a 



+ 



are the usual creation and annihilation op- 



erators [eta), o,7.i] = 2ttS(u> — to'). The A = vacuum state 
is defined by a w |0) = 0. 

For our purposes it will be sufficient to assume that at 
A = the quantum field is in its vacuum state [pi. To 
study the time dependence of the theory it is best to work 
in the Heisenberg representation where the operators are 
functions of time and the quantum state does not change. 
The time evolution of the Heisenberg fields are given by 
the wave equation Eq. (^) supplemented by the condition 
that 4>q and ttq ~ d(f>o(r])/dX reduce to their initial val- 
ues on the A = surface. While this is a complicated 
equation it will be sufficient, for our purposes, to analyze 
it in the geometric optics approximation. Eventually all 
formulas for physical quantities will be expressed geo- 
metrically, so they will not depend upon the coordinate 
system. To explain the geometric optics approximation 
however, it is simplest to work in Painleve coordinates; 
i.e., A, r, since the dependence of the solution on A and r 
factorizes in these coordinates. 

If we look for solutions of the form <fio = r _1 e lwA / w (r), 
it is easy to see that, for large u>, f u (r) can be written as 



Uf u (r) =iojS 1 , 2 (r) + 0(uj- 1 ), 
Si, 2 (r) =±r-2Vra^±2aln| 



(9) 



/i,a(aO = 



Ai(0±7Tl(0T 



l(0 



£3/2 



(12) 



(f>[ = d0i/d£ and 51,2(^1,2) = X + Si t2 (r). Now, since the 
field 4>\ and the momentum 7ri are expressed through 
creation and annihilation operators defined at A = 0, the 
above set of equations allows us to compute any Green's 
function of the field O a t any later time. 

A quantity which is of special interest is the expecta- 
tion value of the A, 77 component of the energy-momentum 
tensor, since it appears on the right-hand side of the Ein- 
stein equations. If the integral of this quantity over a 
surface of fixed r approaches a constant as r — > 00, this 
implies a constant flux of outgoing energy, i.e., Hawking 
radiation. We wish to emphasize that the technique de- 
scribed below permits us to calculate T\ v for any finite 
values of A and rj but for the purposes of this Letter we 
will concentrate on the energy flux through an infinitely 
large sphere. 

In Lemaitre coordinates the total energy flux through 
an infinitely large sphere is given by 



J 



4tt 



r 5/2 
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Xri) 



(13) 

where (Txrj) is the vacuum expectation value of the off- 
diagonal component of the stress-energy tensor of the 
field 4>o ( a normalization factor of (47r) _1 is introduced 
since <po denotes the 5*-wave component of the massless 
field) . Using the explicit expression for the time evolution 
of the field 0q , Eq. (|l^) , it is a straightforward matter to 



3 



compute this flux. In order to do it we need to regularize 
intermediate expressions and we choose to regulate T\ v 
by point-splitting. To this end we introduce a — X 2 — Ai 
and define 



(T\ n ) = lim 



1 



2 



^o(Ai,ry) d(f> Q (X 2 ,r]) 



dX 



dr] 



10), (14) 



where {A, B} denotes the anti-commutator of the fields. 

The calculation of this quantity for large values of A 
and n is simplest in Painleve coordinates. Since the value 
of 77 is fixed, the equation A2 — Ai = a implies that n 
and r 2 are related by 



r-2 = n 



.1/2 



aa 
4rf 



lA3 



+ 0(a 4 ). (15) 



Taking the derivatives in Eq.(|lJ) and considering the 
limit ri t 2 — ► 00, we derive the following expression for 
the flux 



J = -\ lim l [<{/( {x x ),f[(yi)})Wi(r 2 , ri , Xl , Vl ) 



+({f2(x2),f 2 (y2)})W 2 (r 2 ,r 1 ,x 2 ,y 2 )}, 
where the functions fi j2 are defined in Eq.ffil 



(16) 



Wi(r 2 ,n,x,y) 



1 



Sl(x)Sl(y) v vn 

and xi,2, 2/1.2 arc defined through the set of equations: 

X 1 +S 1 (n)=S 1 (x 1 ), X 2 + S 1 (r 2 ) = S 1 (y 1 ), (17) 

and the same for x 2 , y 2 with the change S\ — > S 2 . Given 
these equations it is easy to derive the relation between 
yx )2 and x\ 2 as a power series expansion in a. It is clear 
from Eq.@ that the points Xi,x 2 are those points on 
the A = surface from which the null-geodesics which 
wind up at the point A, 77 must start out. From the spe- 
cific form of the functions S±^ 2 we see that in the limit of 
large A and large 77, the limiting values for these points 
are x\ — > 00, x 2 /a — > 1. 

Let us now compute the anticommutator 
{{fi(xi),fi(yi)}). Starting from Eq.@ we derive 



({/i(zi),/i(yi)}>= 

+({4>' 1 (x 1 ),4>' 1 (y 1 )})- 



(18) 



where the dots represent the terms which do not con- 
tribute to the flux in r± — ► 00 limit. The anticommuta- 
tors in the above formula are computed using Eq. (^) and 
we find the energy flux 



({7Ti(^i),7ri(yi)}) = 



-2 (X. 



7T 



x\r 



4xij/i 



?r(xf - x\) 



2\ 2 ■ 



(19) 



so that the final result for the anticommutator reads 

1 1 



({f[(x 1 )J[(y 1 )}) = 



2n (xi ~ yi)'- 



+ 



(20) 



Performing a similar calculation for the second term in 
Eq.(|l6|) and substituting expansions of 2/1.2 and r 2 in 
terms of xi 2 and n, we find 



J 



1927ra 2 12 



(21) 



where we have introduced Hawking temperature Tjj = 
1/(8ttGM). 

We should make two comments concerning this result. 
First, the energy flux at large distances, we have just de- 
rived, is finite. This is in accord with the expectations 
regarding possible divergences in the stress-energy tensor 
computed in a gravitational background || . Quite gener- 
ally, the result of such a calculation is divergent and the 
divergences are usually removed by appropriate renor- 
malization. The allowed counterterms have a restricted 
form and include renormalization of the cosmological and 
Newtonian constants, as well as some other terms that 
do not appear in Einstein's action. It turns out that none 
of allowed counterterms can renormalize the off- diagonal 
element of the energy momentum tensor in the gravita- 
tional background of a Schwarzschild black hole, both in 
Schwarzschild and Lemaitre coordinates. Therefore the 
result for the energy flux should come out finite and, as 
we have seen, it does. The second comment concerns the 
calculation of the energy flux through the sphere of finite 
radius. It turns out that the result we obtain in that 
case is not finite in that lnu terms remain. We believe 
this to be due to the fact that we have restricted atten- 
tion to spherically symmetric configurations and have not 
considered higher angular momenta. This is supported 
by the fact that in the case of a two-dimensional black 
hole, where higher angular momentum modes are absent, 
our result for the flux is finite for arbitrary values of r. 
We should emphasize that both the outgoing and the in- 
falling fluxes are separately divergent and only the sum 
of the two provides a finite, unambiguous answer. 

To see that this time-independent flux of energy at 
large distances corresponds to what one would expect 
from a body at a Hawking temperature Th it is nec- 
essary to weakly couple the massless field to a detector 
|^|, 0] (which acts as a thermometer) located at some fixed 
Schwarzschild radius r. In order to make the computa- 
tion more realistic we consider adiabatically switching on 
the detector at some time t = to for a finite amount of 
time 5. (Note, adiabaticity requires E typ S ^> 1, where 
E typ is typical detector level spacing.) To realize this 
situation we add an interaction term to the free field 
Lagrangian of the form V int ~ e -(*-to) 2 /(2<5 2 ) ( f )0 (t, r )M, 

where M is an operator which acts in the Hilbert space 
of detector eigenstates. It follows from second order per- 
turbation theory that the probability of exciting the de- 
tector to a state of energy E is proportional to the Fourier 
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transform of the Green's function of the massless field || : 

V(AE) ~ \(E\M(0)\E Q )\ 2 J &t&t' e~ lAE{t -^ x 

e -[(t-to)>+(f-to)»]/(M a ) <0o( t )r )^ (t ' jr ) >j (22) 

where AE = E — Eo and Eo is the ground state energy 
of the detector. As in the calculation of the flux which 
we already described, the Green's function in Eq. (|22|) is 
computed using the evolution equation for the field </>o 
which relates it to initial conditions for (f>o and ttq on the 
surface A = 0. It is convenient to define £1,2 and 
as in Eq.(|l7|) and identify n = r-i = r, Ai = Xi(t,r) 
and A2 = \i{t' ,r). As before, these are the points from 
which infalling and outgoing geodesies leave the A = 
surface in order to reach the points t,r and t',r. Now 
consider Eq. (|2^) for large values of t . Using the explicit 
form of the functions Si 2 it is easy to find the following 
approximate solutions: xi = t + r, yi = t' + r, X2 = 
a(l + 2 e -(*-'-)/( 2Q )), y 2 = a (l + 2 e -(*'- r )/( 2Q )), with 
t ~ t' ~ to- It is then clear that, asymptotically, X2 — > 
t/2 — ► a and x\ — > y\ — > 00. In this limit, the Green's 
function can be written as 



-1 



47rr 2 



(In |xi - 2/1 1 +ln \x 2 - y%\ 



+ y [k(xi,vi) + n(y 2 , x 2 )] + c ) , 



(23) 



where ft(a;, y) = 9{x — y) — 9{y—x) and c is some constant. 

It is instructive to consider the terms in Eq.(^||) sepa- 
rately. The constant does not contribute to V{E) since 
it yields a result proportional to exp(— AE 2 8 2 ) <C 1. The 
In \xi — y\ I term and the terms described by the function 
k give simple, a independent, contributions that can be 
written as 



Vi(AE) 



a~e 



■OdSAE)- 1 ). 



(24) 



The important part of the final result comes from the sec- 
ond term in Eq.(p3|) which describes the radiation coming 
from the vicinity of the horizon. Appropriately shifting 
the integration variables we obtain 

V 2 (AE) ~ - J dtdt ' e -^E(t't') e --[(t-t r+(t'-to) 2 ]/(2S-) 

xlnle-'/^-e-*'/^)!. (25) 

If we then change the variables to v = (t+t 1 ), u = (t—t 1 ) 
and neglect all the suppressed terms we arrive at 



1 



1 



,a_e/t h _ 1 2 



Since the total probability is given by the sum of Vi and 
V2 we get the final result: 



V(AE) I <^7|7V^(0)|^7 > | s 



1 



(26) 



S ' AE - e A£/T H _f- ( 2? ) 

The interpretation of this formula is straightforward. If, 
at a large distance from the black hole, an observer 
switches on a detector which interacts with the mass- 
less field for finite amount of time, then the energy levels 
of the detector get populated as if the detector was in 
equilibrium with a thermal distribution of particles at a 
temperature Tjj. 

In conclusion, in this Letter we have presented a canon- 
ical Hamiltonian derivation of Hawking radiation for the 
case of a Schwarzschild black hole. In this formalism 
the Hamiltonian of the system is perforce time-dependent 
and so the discussion of the vacuum state of the massless 
field is irrelevant to the long time behavior of the system. 
By the very nature of our construction, the state of the 
system is always a pure state and for this reason there 
are correlations of the fields inside and outside the black 
hole horizon. We also explicitly demonstrated, in agree- 
ment with arguments that concern the general structure 
of the divergences in the energy- momentum tensor in the 
gravitational background, that the Hawking flux at infin- 
ity is a uniquely defined finite quantity. This opens up 
the possibility of identifying the time-dependent func- 
tion that should be inserted into the right hand side of 
the Einstein equations to study the back reaction of the 
Hawking radiation on the geometry of the space-time in 
a self-consistent way. 

Clearly, this approach allows us to address more issues 
than we can touch upon in a Letter; in particular, the 
question of back-reaction. Such issues will be discussed 
in detail in a long paper which is in preparation. It would 
be remiss of us, however, to conclude without pointing 
out that there is a significant difference between the case 
of the two and four-dimensional black hole. Direct sub- 
stitution shows that the solutions in Eq. (0) are exact for 
the two-dimensional black hole. This is not true for the 
four-dimensional black hole where the geometric optics 
approximation breaks down at r — > 0. Thus, the be- 
havior of the quantum theory near r — > is different in 
the two cases. This is important because study of both 
cases shows that after a finite time some modes of the 
field no longer appear in the time-dependent Hamilto- 
nian and thus, their time evolution becomes trivial. It 
requires careful study of the behavior of the system in 
order to understand whether or not these modes stay de- 
coupled as the hole evaporates and whether they store 
any information. 
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